Static and Dynamic Correlations in Many-Particle Lyapunov Vectors 
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We introduce static and dynamic correlation functions for the spatial densities of Lyapunov vector 
fluctuations. They enable us to show, for the first time, the existence of hydrodynamic Lyapunov 
modes in chaotic many-particle systems with soft core interactions, which indicates universality of 
this phenomenon. Our investigations for Lennard- Jones fluids yield, in addition to the Lyapunov 
exponent - wave vector dispersion, the collective dynamic excitations of a given Lyapunov vector. 
In the limit of purely translational modes the static and dynamic structure factor are recovered. 
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The dynamical instability of systems, which is quan- 
tified by the spectrum of Lyapunov exponents, has long 
been recognized as a most important characteristic for 
chaotic systems Q. Its understanding for extended dy- 
namical systems provides essentials for the foundations 
of equilibrium and non-equilibrium statistical mechan- 
ics 2]. This explains the on-going scientific activities 
aiming at the connection between Lyapunov spectra of 
many-particle systems and macroscopic properties such 
as transport coefficients Q or entropy production rates 
Q S To each of the Lyapunov exponents in the 
spectrum one associates a Lyapunov vector in state or 
phase space, which in contrast to the Lyapunov expo- 
nent is a time-dependent quantity. Apart from singular 
early investigations Q these mutually orthogonal vec- 
tors were for a long time used mainly as a means for 
calculating Lyapunov spectra. The scientific interest in 
these quantities raised when it was realized that espe- 
cially in extended systems they contain non-trivial infor- 
mation about the perturbation dynamics. For instance, 
for certain coupled map lattices it was found that the 
Lyapunov vectors are often localized in state space @, 
a fact that was later also observed in many-particle sys- 
tems 0,0,^3- Very recently, however, a very surprising 
phenomenon was observed. In contrast to the localized 
behavior of vectors belonging to the largest Lyapunov ex- 
ponents, it was found that the Lyapunov vectors corre- 
spondingto near-zero exponents exhibit wave-like struc- 
tures |rl Ttol ITU . Due to the approximately linear rela- 
tion between Lyapunov exponent and wave number of the 
associated vector the latter were termed hydrodynamic 
Lyapunov modes. These findings immediately triggered 
theoretical efforts to gain deeper insights on the basis of 
simplified models such as products of random matrices 
|T3 |. However, a really satisfactory state of understand- 
ing has not yet been achieved. For instance, molecu- 
lar dynamics simulations indicated that hydrodynamic 
Lyapunov modes seem to exist only in hard-sphere or 
hard-disk fluids, but not in many-particle systems with 
soft interaction potentials [y, flCM. Illj . Although there are 



some arguments for this distinctive behavior based on 
fluctuations of finite time Lyapunov exponents, the ex- 
istence of such a sharp difference between hard and soft 
potentials remains surprising. The purpose of this Letter 
is to show that hydrodynamic Lyapunov modes actually 
do exist also in soft-potential many-particle systems such 
as Lennard- Jones (LJ) fluids. This is achieved by con- 
sidering static and dynamic correlation functions for the 
density of Lyapunov vector fluctuations, which we define 
appropriately in the spirit of generalized hydrodynamics 

m 

Before we present our results for Lennard-Jones fluids 
let us briefly recall the definition of Lyapunov vectors and 
their dynamics. For an autonomous dynamical system 
in M. D given by T — G(T) the tangent space dynamics 
describing small perturbation around a reference trajec- 
tory T(t) is given by ST = M (T(t)) ■ ST with the Jacobian 
M = For a D— dimensional dynamical system there 



exist D Lyapunov vectors e^ a '(t) £ MP, a 
which obey the equations of motion |14L Il5| 



L. 



! (a) = M(T(t)) • e (Q) 
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with m^t) = [e^ ■ M(T(t)) ■ + • M(T(t)) ■ eW] 
x(l — 5 a p/2). Eq.0J describes the dynamics of a frame 
consisting of D perturbation vectors ST which are kept 
orthonormal by the nonlinear terms containing the 
Lagrange multipliers m a p(t). The Lyapunov exponents 
A 1 -"' are given by the time averages A'"* 1 = {m aa (t)) t , 
with (. . .) t = lim t ^oo j J_ . . . dt' . In practice one usually 
does not integrate Eq. Jjj because of numerical instabil- 
ities, but rather obtains the Lyapunov vectors e^ a \t) 
(and Lyapunov exponents) by the so-called standard 
method consisting of repeated Gram-Schmidt orthogo- 
nalization of an evolving set of D perturbation vectors 
ST(t) 1\. One should note that with this procedure the 
Lyapunov vectors are automatically ordered such that 

belongs to the largest Lyapunov exponent 
to the second largest A*- 2 -*, and so on. Eq.Q corresponds 
to the continuum limit of this procedure and is given 
here explicitly to emphasize that e' a ' (t) is a well defined 
dynamical variable in continuous time. 

We treat the Hamiltonian dynamics of N interacting 
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FIG. 1: Snapshot of the LV density u {a) (r,t) for a 1-d 
Lennard- Jones fluid (TV = 100) for which the corresponding 
Lyapunov exponent A^ 96 ^ is close to zero. The displacements 
lif' drawn vertically at the current particle positions Ri show 
indications of a wave-like structure. The insert shows for this 
system the positive branch of the Lyapunov spectrum. 
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FIG. 2: Contour lines of the static correlation functions 
S (aa) (k) for a 1-d Lennard- Jones fluid (N = 100, L = 1000, 
T = 0.2) plotted as function of k and X^ a K The ridge at low 
k— and A— values as indicated by the dashed line, implies the 
existence of hydrodynamic Lyapunov modes. 



particles of a d— dimensional fluid. Therefore T = (R, P) 
is a phase space vector in D = 2dN dimensions evolving 
as r(<) = exp(iCt) T(0) where C denotes the Liouville 
operator. The coordinate part of T is written as R — 
. . . , Rn), with Ri £ R d being the coordinate vector 
of the i—th particle. The momentum part P is defined 
analogously. Correspondingly the Lyapunov vectors 
can also be split into a coordinate and a momentum part 
as e( Q ) = («i , . . . ,Ujif,Vi , . . . ,Vpf). In the follow- 
ing we will concentrate on correlations in the coordinate 
part of these vectors. For this purpose we introduce the 

dynamical variable uW(r, t) — J2f=i u ^{^) ~~ Ri(t))i 
which describes the spatial density of the coordinate part 
of Lyapunov vector fluctuations, which we call for short 
LV density. This variable is defined in analogy to other 
microscopic density fluctuations, e.g. of the particle cur- 
rent, in molecular hydrodynamics |l3| . A snapshot of 
this dynamical variable is shown in Fig.l for a 1-d LJ 
fluid. For translational invariant systems as treated here, 
it is convenient to consider the correlations of its Fourier 
transform (k, t) = Ju^ a \r,i) exp(— ikr) dr given by 

N 

u< a >(M) = $>i a) W exp(-ikRi{t)) (2) 

i=l 

The intermediate two-time correlation function of this 
variable is defined as 

F {aa \k,r) ^N- 1 (u {a \k,t + r) u (q > (-*,*)) (3) 

which for d ^ 1 is a second rank tensor. From the 
latter we get for r = the static LV density corre- 
lation function S( QQ )(fc) = F^ aa \k,T = 0) and by 
Fourier transformation with respect to time the dy- 
namical LV density correlation function S^ aa '(k, u>) — 
(27T)- 1 / F < - aa 1(k,T)exp(iujT) dr. Note that in contrast 
to the Lyapunov vectors itself, their correlation functions, 
being time averages, define global quantities similar to 



the Lyapunov exponents . Even more far reaching is 
the following fact. It has been shown [l^ that generi- 
cally for long times the Lyapunov vectors e^ a \t) depend 
on time only via its current phase space point T(t), espe- 
cially it does not depend on the initial conditions e^(t — 
0) as Eq.Q might suggest. This implies that in long time 
averages e' a '(i) can be replaced by e^ a \T(t)), where the 
vector fields e( Q )(r), a = 1, . . . , D, provide what in 0] 
is called a " stationary Lyapunov basis" . Accordingly the 
evolution of the dynamical variables, e.g. Eq.Q, can be 
expressed as u ( a )(k,t) = exp(iCt) uW(Jfe,0) and, concep- 
tually most important, the time averages in our correla- 
tion functions can be evaluated also as ensemble averages 
F( QQ )(fc,r) = N- 1 (u^(k,t) u( Q )(-fc,0)) over the in- 
variant density /i(r), i.e. (...) = J . . . n(T) dT. Further- 
more standard techniques for the analytical evaluation 
of these correlation functions, such as the Mori-Zwanzig 
projection operator formalism [l3| . could be applied in 
principle. The problem, however, is the lack of knowledge 
about the vector fields e^ Q ^(r). One of the few known 
features is that for systems with translational invariance 
there exist Lyapunov vectors with constant spatial com- 
ponent = c^" ) belonging to the zero Lyapunov 
exponents A^ Q °^ = related to this symmetry n^. In 
these special cases the LV density variable u^ a °'(r,t) is 
proportional to the ordinary microscopic particle density 
«( Q »'(r,t) = c^p(r,t) with p(r,t) = S ( r ~ ^(*)) 
|17|. Correspondingly for these zero modes the compo- 
nents of the LV density correlation matrices S^ a ° a °^{k) 
and 5( aoQ °)(fe , u>) are simply proportional to the static 
and dynamic structure factor S(k) and S(k,ui), respec- 
tively. 

The general behavior of S^ aa \k) and S^ aa \k,u)) is in- 
ferred from numerical simulations. The particle dynam- 
ics is obtained by integrating Hamilton's equations with 
periodic boundary conditions using a standard truncated 
LJ interaction potential given by V(r) — 4e[(<j/r) 12 — 
(a/r) 6 ] — V c for r < r c — 2.5er and V(r) — elsewhere 
|l8| . V c is chosen such that V(r) is continuous at r c , 
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and the parameters a, e and the particle mass m are 
all set to unity. This means that e.g. system length L, 
(kinetic) temperature T and time is measured in units 
of cr, e, and (r7icr 2 /48e) 1 / 2 , respectively. The Lyapunov 
vectors e^ a \t) were obtained iteratively with the stan- 
dard method pj. The time step used in the symplectic 
Verlet algorithm of our MD simulations is A = 0.008 
and the re-orthogonalization for the calculation of the 
e^ a \t) is repeated in time intervals of 30 — 100A. The nu- 
merical evaluation of the correlation functions is started 
only after appropriate equilibration periods. The coor- 
dinate part of a Lyapunov vector, as shown in Fig.l for 
L = 1000 and T = 0.2 is a quantity fluctuating in space 
and time. The equal time correlations in space are cap- 
tured by S( aa \k). For simplicity we present first results 
for the one-dimensional LJ fluid, for which S' aa ' (k) is a 
scalar quantity. In this case we can depict its behavior 
in a contour plot. Since for A'"- 1 > the latter decreases 
monotonically with the index a, we can plot S^ aa \k) 
also over the A — k— plane thereby introducing the func- 
tion S(X,k). The surmise, which is supported by our 
numerical simulations, is that the latter converges with 
increasing N (for fixed temperature T and density N/L) 
to a well-defined limit. In the contour plot of Fig. 2 we 
recognize a ridge at small A and k, which defines a dis- 
persion relation A(fc) or fc(A) implying that the Lyapunov 
mode corresponding to Lyapunov exponent A is domi- 
nated by a spatial oscillation with wave number k(X). 
In Fig. 3 we determined the dispersion relation A(fc) by 
plotting A versus fc max = argmax^ S(X, k) for various 
systems differing in density and temperature. We find 
that the dispersion X{k) appears to be independent of 
these parameters in the regions where X{k) is well de- 
fined. Only the extent of the latter regions is parame- 
ter dependent, as it decreases e.g. with decreasing den- 
sity. The common dispersion relation X(k) being linear 
on a log-log scale results in a power law X(k) ~ k v with 
T] = 1.2 ± 0.1 although a linear dispersion with quadratic 
corrections cannot be excluded. In any case these re- 
sults unambiguously show the existence of hydrodynamic 
Lyapunov modes in LJ fluids. Actually, our investiga- 
tions show that their existence is not restricted to 1-d 
systems. In isotropic fluids with d > 1 the Cartesian 
components F^ o) (jfe,r) of F^ aa \k,T) can be written in 
terms of longitudinal and transverse correlation functions 
F ( L aa) and 4 aa) as F^ a \k,r) = k^KF^Qk] ,r) + 
(<W - knk)F^ aa \\k\,T) with k„ = (k/ |fc|V We find 

for d = 2 that F ( L aa \\k\ , 0) and F<f a) (\k\ , 0) behave sim- 
ilar to the 1-d case shown in Fig. 2 implying the existence 
of hydrodynamic Lyapunov modes also there. 

More detailed information can be obtained from the 
dynamical LV correlation functions S^ aa ^ (k, w), which 
encode in addition to the structural also the tempo- 
ral correlations. As usual the equal time correlations 
can be recovered by a frequency integration S^ aa \k) = 
j S^ aa \k,u>) dm. In Fig. 4 we show typical examples 
for S^ aa '{k,w). It consist of a central "quasi-elastic" 
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FIG. 3: The dispersion relation X(k) of the hydrodynamic 
Lyapunov modes shown for different particle densities and 
temperatures, appears to be independent of the latter for a 
wide range of parameters. 
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FIG. 4: The dynamic correlation function S^ aa ^ (k,aj) for a = 
96 and k = 2n/L. The full line results from a fit by a 3-pole 
approximation. The corresponding decomposition into three 
Lorentzians is also shown. The insert shows a comparison of 
such fits for k = 2n/L, 4-7r/L, and Stt/L (from the top). 



peak with shoulders resulting from dynamical excitations 
quite similar to the dynamical structure factor S(k, u>) of 
fluids. In order to extract the dynamical information 
we used a 3-pole approximation for the Laplace trans- 
form of i^( aa ) (fc, r) from Eq.©, which amounts to fitting 
S^ aaS> (k,uj) by a superposition of three Lorentzians, one 
central peak at lo = and two symmetric peaks located at 
uj = ±Lu(k). The fits are also shown in the figure and de- 
scribe the frequency dependence of S^ 010 ^ (fc, lo) quite well. 
Such a dependence arises naturally e.g. from continued 
fraction expansions based on Mori-Zwanzig projection 
techniques [l3j, which as mentioned may be applied also 
to this problem. These fits allow us to extract the disper- 
sion relations a/"- 1 (k) for each of the hydrodynamic Lya- 
punov modes with index a. The result is shown in Fig. 5 
for several of the Lyapunov modes. Clearly, this tells 
us that a Lyapunov mode corresponding to exponent A 
is characterized, apart from the dominating wave length 
k(X), by a typical frequency w(fe(A)). Because ^ is non- 
vanishing, this implies propagating wave-like excitations. 
The origin of the characteristic frequency u>(k(X)) is not 
yet fully understood. Probably it reflects rotational mo- 
tion of the orthogonal frame {e' a )(£)} around its refer- 
ence trajectory |6|. The full LV dynamics, however, is 
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FIG. 5: Dispersion relations uj^ a '(k) (top) and the 
k— dependence of the width of the central peak (bottom) for 
different Lyapunov vectors obtained from 3-pole approxima- 
tions as in Fig. 4. Note the clear separation of the correspond- 
ing time scales. 

more complex. We find, for instance, that the coherent 
wave-like motion is switched on and off intermittently 
|19|| . This phenomenon presumably is responsible for the 
strong central peak in S^ aa \k,u) (see Fig. 4). 

In summary, we have shown that the LV density cor- 
relation functions provide a valuable tool for quantifying 
spatiotemporal phenomena in many-particle Lyapunov 



vectors. They allowed us to identify for the first time 
hydrodynamic Lyapunov modes and their dynamics in 
soft-potential systems. These findings indicate that the 
existence of hydrodynamic Lyapunov modes is a univer- 
sal phenomenon in chaotic translation invariant many- 
particle systems, independent of the interaction poten- 
tial. This view is supported by the fact that these modes 
appear to be present also in one- and two-dimensional flu- 
ids with Weeks-Chandler-Anderson (WCA) interaction 
potentials |2jj. There exist obvious generalizations of 
our correlation functions, which will be treated in forth- 
coming papers. For instance, one can consider corre- 
lations between different Lyapunov vectors resulting in 
cross-spectra S^ a ^'(k,iv), or correlations of the LV den- 
sity not only in coordinate space but in phase space 
e^(r,p,t) = Elicit) Sir-R^Sip-P^t)), which 
is an approach used similarly in kinetic theory [Lj • 
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